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Abstract 

The virtues of an effective field theory (EFT) approach to many-body prob- 
lems are illustrated by deriving the expansion for the energy of an homoge- 
neous, interacting Fermi gas at low density and zero temperature. A renormal- 
ization scheme based on dimensional regularization with minimal subtraction 
leads to a more transparent power-counting procedure and diagrammatic ex- 
pansion than conventional many-body approaches. Coefficients of terms in 
the expansion with logarithms of the Fermi momentum are determined by 
the renormalization properties of the EFT that describes few-body scatter- 
ing. Lessons for an EFT treatment of nuclear matter are discussed. 
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I. INTRODUCTION 



The effective field theory (EFT) approach exploits the separation of scales in physical 
systems The application of EFT methods to many-body problems promises insight 

into the analytic structure of observables and a consistent organization of many-body cor- 
rections, with reliable error estimates. Our ultimate goal is to calculate the properties of 
nuclear matter and finite nuclei, but there are many issues as yet unresolved concerning 
the practical application of EFT methods to these systems. As a step toward this goal, we 
consider a simpler problem for which we can illustrate some essential features that might be 
obscured in the full nuclear matter treatment. 

The dilute Fermi gas has often served as a benchmark many-body problem, despite a 
shortage of physical applications. Recently, there has been experimental progress in creating 
and measuring dilute low-temperature Fermi systems 0, although these experiments at 
present are not sensitive to interactions. In this paper, we consider a uniform system of 
nonrelativistic fermions of mass M, interacting via a purely repulsive, spin-independent 
underlying interaction with characteristic range R = 1/A. The restriction to repulsive 
interactions allows us to postpone discussion of Cooper pair and liquid-gas instabilities 

The prototype interaction of this type is a hard-sphere repulsion of radius R. For mo- 
menta p <^ A = 1/R, two-body scattering with a hard-sphere potential is described by the 
effective range expansion [see Eq. (j^)]. All of the effective range parameters are "natural"; 
that is, their magnitude is given by a power of the underlying scale R times a dimensionless 
number of order unity. For general interactions, natural values are expected for the effective 
range and higher-order parameters, but not necessarily for the scattering length, which can 
be arbitrarily large if there is a bound state near zero energy. The latter case is of particular 
interest, of course, because of the nature of s-wave nucleon-nucleon scattering. However, a 
fine-tuned scattering length requires special treatment so our discussion here will 

assume natural scattering lengths. 

To our knowledge, the most complete low-density expansion for the ground-state energy 
per particle of such a system is |]T2|-|I^: 
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In Eq. (|l]), tts and are the s-wave scattering length and effective range, and Op is the 
p-wave scattering length. The spin degeneracy is denoted by g. For a natural system, this is 
an expansion in Fermi momentum kp over the scale A. The mean-field correction of 0{kp) 



dates from 1929 [15|, the 0{kp) correction from the 1950's |16|, while the 0{kp) corrections 



and the logarithm were found in the 1960's fl^. The complete expression in Eq. (|I|) has been 



derived using the method of correlation functions |T2|, by expanding Goldstone diagrams 
T3| , p^ , and by expanding Feynman diagrams fl^. Here we rederive and illuminate this 



result using EFT methods. 



While field theory approaches to many-body physics are well-established {e.g., Refs. |T7 



and [|l^]), the treatment of divergences and renormalization issues for a dilute Fermi system 



I^^14[ are reminiscent of the old-style approach to renormalization theory, in which renor- 



malization is primarily a technical device to get rid of divergences in perturbation theory 



{i.e., to sweep problems under the rug) |[T9|. The modern viewpoint in particle physics that 
"renormalization is an expression of the variation of the structure of physical interactions 
with changes in the scale of the phenomena being probed" |T^ has led to powerful and 
elegant renormalization group methods and to the development of the effective field the- 
ory approach, which provides low-energy descriptions in terms of local, nonrenormalizable 
interactions. 

We want to revisit the dilute Fermi gas with the same type of insight that effective field 
theory has given to particle physics. Some first steps were taken in Ref. in which a 
double expansion in momentum and the strength of a model potential was exploited to 
study renormalization and EFT breakdown at finite density. Here we extend the analysis 
to a general low-density expansion that is perturbative only in the Fermi momentum. 

Effective field theories for scattering of two fermions in the vacuum are now well under- 
stood and are particularly straightforward for systems with natural scattering lengths 
2T|,|22|]. At low energies, with momentum k <^ A, the underlying interactions are not re- 



solved, so they may be replaced by a series of local (contact) interactions; this description is 
completely general to the order of truncation. Contributions to the scattering amplitude are 
expressed as a systematic diagrammatic expansion in k/A. The prescription for determining 
what diagrams contribute to a given order in the expansion is called power counting. Since 
the contact interactions are singular, they require regularization and renormalization, and 
the power counting prescription depends on the subtraction scheme. It is particularly useful 
to have a subtraction scheme in which any Feynman graph with a particular vertex only 
contributes to the expansion of the amplitude at a particular order, and for which power 
counting is simple dimensional analysis. For a natural system, these conditions are satisfied 



by minimal subtraction implemented using dimensional regularization 1 2 1,22 



The EFT expansion for two-body scattering efficiently reproduces the effective range 
expansion, which is a Taylor expansion in external momentum p divided by the characteristic 
scale of the underlying interaction A. When using a low resolution probe, only the scale A of 
new, unresolved physics enters. Up to factors of order unity, the radius of convergence of the 
expansion {i.e., where errors are under control) is A, even though the domain of analyticity 
may be infinite. A reliable estimate for the truncation error is order unity times k/A to the 
first omitted power. 

For a problem like this where the physics is "transparent", the calculation would ide- 
ally be equally transparent. The EFT treatment in the vacuum does exactly that. One 
might expect a simple analysis to be possible at finite density as well, but the standard 
treatments are actually quite complicated, with summations to all orders and expansions 
for each diagram The EFT approach presented here maintains the virtues of the 

EFT description in the vacuum: each Feynman diagram contributes to a definite order in 
the kp/A expansion of the energy. Previous work strained to eliminate divergences; with an 
EFT they are natural, understood, and even exploited. 

The analytic structure of the low-density expansion as a function of kp is dictated by 
the interplay of short- and long-distance behavior. This is particularly well manifested in 
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the effective field theory approach. Recently, Braaten and Nieto showed how renor- 
malization in a dilute boson gas determines the coefficients of logarithms in the low-density 
expansion of the energy. Here we make an analogous analysis for fermion systems to find the 
coefficient of the logarithm in Eq. (|l|) from renormalization group arguments. Extensions 
of this analysis may be fruitful in deriving the analytic structure of the energy functional 
of nuclear matter and other strongly interacting fermion systems. We also discuss how to 
complete the renormalization, which requires additional input beyond what is determined 
by two-body scattering. 

In section II, we review how EFT power counting in the vacuum generates the effective 
range expansion. In section III, the EFT is used to reproduce the analytic terms in Eq. (|1]) 
and in section IV the logarithm is discussed. We focus on the physics with a minimum of 
formalism. More complete details will be given elsewhere, including numerical details and 
extensions to finite temperature and other observables 0. The results are put in the context 
of the nuclear matter problem in Section V. 

II. EFT AT ZERO DENSITY 

In this section, we review the EFT for two-body scattering in the vacuum.Q We consider 
a natural EFT for heavy nonrelativistic fermions of mass M, whose interactions are spin- 
independent and with strength correctly estimated by naive dimensional analysis based on 
an underlying scale A. In the center of momentum frame, the scattering amplitude T of two 
particles with momenta of magnitude k and scattering angle 6 can be expanded in partial 
waves as 

where Si is the scattering phase shift for angular momentum / and Pi is a Legendre poly- 
nomial. Furthermore, for short-range interactions, k cot 6i can be expanded in powers of /c^ 
as 

11 3 
kcot5o = \- -r^k'^ + . . . , kcot5i = -—^ + . . . , (3) 

where a^, a^, and are the s- and p-wave scattering lengths and s-wave effective range, 
respectively. Inserting Eq. (|) into Eq. (||) and expanding the result up through OikP' / Is?'), 
we obtain 

T(fc, cos d) = [l - ia^k + (a,r,/2 - a^) A;'] - cos 9 + Oik^A^) . (4) 



M 



M 



For a natural system, ~ Op ~ ~ 1/A = i? and we have a perturbative expansion 
in k/A. An example commonly studied is a hard-sphere gas with radius R, in which case 
ttg = ap = R and = 2R/3. 



^For a more complete treatment, see Refs. |21,22| 
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Since we assume <^ A, all interactions in the EFT are short-ranged and we have only 
contact interactions. Thus to reproduce this expansion with an EFT we consider a general 
local Lagrangian for a nonrelativistic fermion field that is invariant under Galilean, parity, 
and time-reversal transformations:^ 



idt 



2M 



2 io L 



+ ^(VW)^ ■ (V'W) - ^(V^V)' + • • • , (5) 

where V is the Galilean invariant derivative and h.c. denotes the hermitian con- 

jugate. The terms proportional to C2 and C'2 contribute to s-wave and p-wave scattering, 
respectively. The dots represent terms with more derivatives and/or more fields. Higher- 
order terms with time derivatives are omitted, as they can be eliminated in favor of terms 
with spatial derivatives pT| , p2| . Relativistic kinematic corrections can be included systemat- 
ically, but all such terms are suppressed by factors of (A/M)^, which we take to be negligible; 
they are considered explicitly in Ref. |p. 

The scattering amplitude for fermions in the vacuum is simply related to a sum of 
Feynman graphs computed according to this Lagrangian. The terms in Eq. involving 
only four fermion fields reduce in momentum space to simple polynomial vertices that are 
equivalent to a momentum expansion of an EFT potential for particle-particle scattering, 

(k'l V^FTlk) = Co + C2(k'' + k2)/2 + k' ■ k + . . . , (6) 

(see Fig. |T| for the Feynman rules). Because of Galilean invariance, the interaction depends 
only on the relative momenta k and k' of the incoming and outgoing particles. The co- 
efficients Co, C2, and C2 can be obtained from matching the EFT to a more fundamental 
theory or to (at least) three independent pieces of experimental data. It is important to note 
that Eq. (|^) is not simply the term-by-term momentum-space expansion of an underlying 
potential because the coefficients also contain short-distance contributions from loop graphs 
(see Ref. pO[ for a simple example of the matching and renormalization of the coefficients). 



The loop integrals that appear in Feynman graphs for two-body scattering with this EFT 
in D spatial dimensions reduce to 

T = f (7\ 

" " i (27r)« A;2 _ ^2 + ■ y') 

For D = 3 spatial dimensions, these integrals are divergent and must be regulated. If we 
apply a sharp momentum cutoff Ac on g, evaluate the integral, and expand in A;, we find the 
leading terms have positive powers of Ac. These are called power divergences. For example, 
when n = 0, 



^The extension to spin- and/or isospin-dependent interactions is straightforward and is illustrated 
in Ref. M. 
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FIG. 1. Feynman rules for (k'|VEFT|k) and the leading 3-body contact interaction. The spin 
indices have been suppressed. 

Lo ^ -^Ac ~^k + 0{k'/K) . (8) 

The power divergences come from the short-distance region of the integral where the loop 
momentum is of order Ac. This contribution is incorrect and depends on the details of how 
we regulate the integral, but can be removed systematically using local counterterms. In 
contrast, contributions from loop momenta of order of the small external momentum k [the 
second term in Eq. (||)] are physical. These contributions are completely decoupled; the Ac 
coefficient tells us nothing about the k coefficient. 

If we imagine two-body scattering mD = 2 dimensions, the situation is different, because 
there will always be a logarithmic divergence (as well as power divergences in general). For 
example, when n = 0, 

Lo^^\n{k/K) + O{l) . (9) 
Ztx 

The dependence on In(Ac) is canceled by a counterterm, but now the contributions from 
the different momentum regions are coupled. Since the logarithms of k and Ac must match, 
the coefficient of the logarithm is independent of the regularization. Further, it can be 
determined by looking only at the log divergence in the integral, without a full calculation. 
We will take advantage of this simplification later when we consider three-body scattering, 
in which logarithmic divergences appear for D = 3. 

Since the power divergences do not provide useful information here and complicate the 
power counting, it is appropriate to use a minimal subtraction scheme, which by construction 
subtracts the power divergences. Dimensional regularization with minimal subtraction is 
particularly convenient for a natural EFT, since loop integrals for two-body scattering pick 
up only the scale of the external momentum (logarithms are also possible when more than 
two fermions scatter). In particular, pT| , p2 
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Thus, the underlying physical scale A appears only in the coefficients and the loops contribute 
only factors of k, which can be determined by dimensional analysis. 

Therefore the estimated contribution of a given diagram is found by counting (i) for every 
propagator a factor M/k"^, (ii) for every loop integration a factor of k^/M, and (iii) for every 
n-body vertex with 2i derivatives a factor fc^YMA^*^^"^^. This scaling of the coefficients 
with M follows from Galilean invariance and with A from dimensional analysis |]2T| , ^ . In 
particular, the vertices with four [n = 2) and six (n = 3) fields scale as 

^''^W^ D,,r.-—-, (11) 

which means that more complicated graphs are suppressed at low energies by additional 
powers of k/A. As above, A is an ultraviolet scale which is determined by the physics 
not explicitly included in the EFT (such as the mass of a heavy exchanged particle). For 
momenta ~ A, the expansion in powers of k/A breaks down and the new physics at the 
scale A has to be explicitly included in the EFT. Therefore, A is also called the breakdown 
scale of the EFT. 

We can use the topological properties of a graph to determine that a diagram with L 
loops or E external lines and Vg" n-body vertices with 2i derivatives scales precisely as k'^, 
where 

oo CO 

n=2 i=0 
n oo oo 

= 5--E + ^^(2z + 3n-5)^"- (12) 

^ n=2 1=0 

At every order in fc/A there are a finite number of graphs that contribute. The mass M 
enters all Feynman diagrams in the same way as a factor 1/M and is irrelevant for the 
relative scaling of diagrams. 

Now we can determine the coefficients by matching. In our case, the experimental 
data we match to are the leading coefficients of the effective range expansion for the scat- 
tering amplitude T, Eq. (Q). Evaluating T in the center-of- momentum frame of the two 
particles, we obtain (see Fig. 0) 



T{k, cos 6) = -Co 1 - i—Cok + -fk' - —Co k' + -^k' cos 6 + 0{kyA') , (13) 
L An Co \47r / Go 

where |k| = |k'| = k and 6 is the scattering angle. Matching Eq. (|13D to Eq. (^), we can 
express the C2i in terms of the effective range parameters: 

Co = ^, C, = Co^, and C', = ^. (14) 
Thus the EFT reproduces the effective range expansion in the vacuum.^ With the interaction 



^If the scattering length is unnatural, the power counting has to be modified and the perturbation 
series must be resummed to avoid a premature breakdown of the EFT. This case has received 
considerable attention recently @-^]. 
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FIG. 2. Expansion of T{k, cosO) in the EFT. 

parameters in the Lagrangian expressed in terms of the effective range parameters, we are 
now ready to apply the EFT to dilute Fermi systems at finite density. 

III. DILUTE FERMI SYSTEMS 



We first review the noninteracting Fermi gas at zero temperature to set our notation 
and conventions. In this case, all states up to the Fermi momentum k-p are occupied while 
the ones above fcp are empty. The density p = N/V of a noninteracting system with Fermi 
momentum kp and degeneracy g = (2s + 1) is obtained by summing over all occupied states. 



P = 9 



d^k 



9{kF - k) 



gj4 
67r^ 



(15) 



In a uniform system, this relationship between p and fcp is unchanged by interactions []17 . 
The energy density Sq of a noninteracting Fermi gas is given by weighting the integral in 
Eq. (|15D with the free single particle energy cj^ = k^/(2M), which leads to 



^0 = Pt 



3 /up 



5 2M 



(16) 



The energy per particle Eq/N is simply So/p. 

Equation (|16D receives corrections from the interaction between the fermions. For low- 
density systems with short-range interactions, p oc <^ l/i?^ ~ A^, and the bulk properties 
of the system are amenable to an expansion in kpR or kp/A. As discussed earlier, the 
EFT method is ideally suited to deal with such a separation of scales. In this section, we 
present an EFT for low-density Fermi systems and illustrate its application by calculating 
the corrections to Eq. ([TgD to order (kpR)^. The use of contact interactions and EFT power 
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counting makes the EFT treatment straightforward and more transparent than traditional 
approaches {e.g., see Ref. [0]). 

The ground-state energy density S can be expressed as a sum of Hugenholtz diagrams, 
which are closed, connected Feynman diagrams with appropriate symmetry factors (see 
Fig. To apply the EFT from the previous section at finite density, we have to adapt the 
power-counting rules from the vacuum to the medium. This is immediate with dimensional 
regularization and minimal subtraction, because the typical momentum flowing in both par- 
ticle and hole lines in energy diagrams is of order kp (see below). The result is clear for hole 
lines, but particle lines can in principle have momenta that are unbounded or as large as an 
imposed cutoff A^. However, as discussed above, the contributions of states with momenta 
close to the cutoff, which depends on the regularization scheme in an EFT, are removed by 
minimal subtraction. The correct effect of high-momentum modes on low-energy physics is 
incorporated in the coefficients of the local operators in the Lagrangian instead. The end 
result is that we can use Feynman diagrams for the ground-state energy, in which holes and 
particles are propagated together, and apply the vacuum power counting with kp playing 
the role of the external momentum. 

Thus we estimate the contribution of a given Hugenholtz diagram by counting (i) for 
every propagator a factor M/kp, (ii) for every loop integration a factor of kp/M, and (iii) 
for every n-body vertex with 2i derivatives a factor fcp /MA^*"*"^""^. Again, the mass M 
enters all Feynman diagrams in the same way as a factor 1/M. Furthermore, the range R 
always has to match up with kp to produce the dimensionless expansion parameter kpR. 
The correct power of R can be inferred from the overall dimension of the diagram, which is 
an energy density. Consequently, it is sufficient to keep track of the powers of kp to classify 
the relative size of diagrams. The power counting rules are conveniently summarized in one 
equation in terms of the vertices of a given diagram. In particular, a graph with V.^ n-body 
vertices with 2i derivatives scales as kp, where u is given by 

oo oo 

i^ = 5 + J2T.i^^ + ^n-5M, (17) 

n=2 i=0 

[this follows from Eq. (|12D with E = 0]. As at zero density, diagrams with more vertices and 
more derivatives are suppressed, and therefore the EFT expansion is efficient and systematic. 

The leading order Hugenholtz diagrams contributing to the energy density are shown in 
Fig. 1^. From Eq. (0), the energy density of the noninteracting Fermi gas is of 0{kp); the 
first correction Fig. ^(a) is of 0{kp). The Feynman rules for the contributions to the energy 
density from the diagrams in Fig. ^ are essentially the same as in standard many-body 
treatments {e.g., see Ref. ||T7|,18|), which we summarize here. 



First, draw all distinct, fully connected diagrams contributing to a given order in kp [as 
determined by Eq. {^IT\)]. Distinct diagrams are those that cannot be deformed to coincide 
with each other, including the direction of arrows. To evaluate a diagram: 

1. Assign nonrelativistic four-momenta (frequency and three-momentum) to all lines and 
enforce four-momentum conservation at each vertex. 

2. For each vertex include a factor —i times the corresponding term from the effective 
potential (k'|VEFT|k), as shown in Fig. |l|. (Recall that k and k' are relative momenta.) 
The spin structure of the vertices has been suppressed in Fig. |l|. For spin-independent 
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FIG. 3. Hugenholtz diagrams for the energy density in a homogeneous, dilute Fermi system. 
The vertices are defined in Fig. 0. 

interactions, the two-body vertices have the structure (Sa-yS/ss + Sas^/s-f), where are 
the spin indices of the incoming hues and 7, 6 are the spin indices of the outgoing hues. 
For spin- dependent interactions one inserts the appropriate operator in spin space at 
each vertex instead. 

For each internal line include a factor iGo{k)a'y, where k = {ko,\<.) is the four- 



momentum assigned to the line, a and 7 are spin indices, and 

' eik-kp) 0{kF-k) \ 



tGoik) 



07 



i6, 



^0 - ^^k + 



(18) 



3. Perform the spin summations in the diagram. In every closed fermion loop, substitute 
-g for each 6 

4. Integrate over all independent momenta with a factor / d^k/{27i)^ where d^k = dko d^k. 
If the spatial integrals are divergent, they are defined in D spatial dimensions and 
renormalized using minimal subtraction as discussed in Sect. || and Ref. |21,22|. For 
lines ending and originating at the same vertex, multiply by exp(ifco'7) and take the 
limit ?7 ^ 0+ after the contour integrals have been carried out. This procedure auto- 
matically takes into account that such lines must be hole lines. 
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5. Multiply by a symmetry factor V('S'ni™(^0'") where S is the number of vertex per- 
mutations that transform the diagram into itself, and m is the number of equivalent 
Z-tuples of lines. Equivalent lines are lines that begin and end at the same vertices 
with the same direction of arrows. 

We illustrate the use of these Feynman rules by applying them to some of the diagrams 
in Fig. |[ We start with the diagram at 0{kp), Fig. |^(a). This diagram has one vertex 
and one equivalent pair of lines originating and ending at this vertex. The overall factor is 
therefore i/2 and we have 

£l = l{-^Co)g{g-l)(^\imJ0^^e^'^^^GoCk)j , (19) 

The dko integration is performed using contour integration, which picks up the 6{kF — k) 
pole, and the remaining d^k integral up to kp is trivial. The result is 

where Eqs. (|Tip and have been used |]T3| . Si is of order kp as expected from the power 
counting. We label the energy contributions by the power of kp that accompanies the generic 
scale R] this power is indicated by the subscript j on £j. 

At Oi^kp) we have two diagrams. Figure ^(c) vanishes at T = because it requires 
particle-hole pairs of equal momentum at the Fermi surface, which do not exist. (In par- 
ticular, there are integrals of the form Jdkf{k)6{kp — k)6{k — kp), which vanish if f{k) is 
regular.) Such diagrams are called "anomalous." Figure ^(b) has S = 2 and two equivalent 
pairs of lines, so the symmetry factor is i/8. Applying the Feynman rules we obtain 

^2 = ^ - 1) / / Goip) Go(p - q) Go{k + q) GS) ■ (21) 

The frequency integrals are again simple contour integrals. A transformation to dimension- 
less center-of-mass variables s, t, and u using 

p = kp{s + t), k = A;F(s — t), and q=kp{t — u) (22) 

leads to 

^2 = p(^7-l)^(A:Fa.)^(0^ Jd's J dH J d'u 6(1 - \s + t\)e{l - \s - t\) 

xO{\s + u\- m\s - u| - 1)^^3^^ , (23) 

where Eq. (0) has been used again. 

The integral over u in Eq. (^3]) contains a linearly divergent term in D = 3 dimensions 
(putting back the factors of kp): 

d^u r°° d^u f^v d^u 

(24) 



fcp t"^ — + ie Jo t"^ — + ie Jo t"^ — + ie 
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When dimensionally regularized, the first term on the right-hand side is purely imaginary 
in our minimal subtraction prescription [cf. Eq. (|l3)]- Since the energy density is real, 
all imaginary parts cancel in the end and the first term does not contribute. We have 
verified this cancellation explicitly. In the conventional approach, one avoids such power- 
law divergences by eliminating the potential, which is implicitly assumed to be valid for 
arbitrarily high momenta, in favor of the scattering amplitude. This amounts to summing 
the "ladder" diagrams to all orders and expanding the result afterwards. The EFT deals 
with such power law divergences more directly by avoiding the assumption of an underlying 
potential. Power law divergences are simply subtracted out and the effect of high energy 
modes is captured in the coefficients of the terms in the effective potential. Renormalization 
theory assures us that this procedure is consistent and correct [||]. As in the conventional 
approach, one could also resum all bubbles with Cq contact interactions in the EFT and 
apply this two-body scattering amplitude instead of the Cq interaction at finite density. 
Since the resummed amplitude falls like 1/p at large momenta, some ultraviolet divergences 
at finite density would be avoided. In the case of nuclear matter, where some mass scales 
are fine tuned to be unnaturally large, such a resummation will have to be done. (See Refs. 
10| , p!I[| for a discussion of this issue in the vacuum case.) For natural scales, however, it is 



much easier to treat all interactions perturbatively. As long as the renormalization is carried 
out consistenly, the results are insensitive to the regularization scheme. 

The second term in Eq. ([2^) is finite and the remaining integrals can be evaluated 



analytically after some judicious partial integrations (see Appendix C of Ref. ||2^ for details). 
The net result for the second-order contribution to the energy density is 

^2 = p(^7-l)^(W^(ll-21n2) , (25) 



which was first obtained by Huang and Yang [|16 . 

At 0{kp), we have seven diagrams contributing to the energy density. Their evalua- 
tion from the Feynman rules is straightforward although somewhat tedious; we will not 
go through the calculations in detail here but only make some comments. Three of the 
diagrams. Figs. H(d), 0(e), and 0(f), vanish when evaluated. This is immediately clear for 
Fig. 0(d), which vanishes due to the same argument as the similar diagram at (9(fcp). In 
Fig. 0(e), the contributions of hole and particle propagation in the line with the tadpole 
loop attached cancel exactly. This cancellation is manifest once the frequency integrals have 
been carried out. Figure 0(f) can again be seen to vanish without any calculation. As- 
sign a four-momentum k to the center loop to which the three tadpole loops are attached 
and perform the dk^ integration. We have to evaluate JdkoG^kY , but only simple poles 
in ko contribute to the integral. The terms with simple poles, however, involve a factor 
6{k — k-p)6{kp — k) times a regular function and so the diagram vanishes. Note that because 
the interactions are contact terms, a large number of exchange diagrams generally vanish as 
anomalous diagrams. 

Diagrams 0(g) and (h), each with three Co vertices, can be directly calculated. The first 
one contains a linear divergence similar to the diagram at (9(fcp), which can be isolated 
{e.g., by transforming to center-of-mass variables), and then removed by minimal subtrac- 
tion, leaving a finite remainder. The second diagram is finite. We have not obtained analytic 
results for these diagrams. After the principal value integrals are carried out analytically. 
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however, the numerical evaluation of the remaining integrals is straightforward. We have 



calculated these integrals using a standard Monte Carlo integration routine [g4| and repro- 
duced the values given in the literature [p!^-p^. Further details will be given elsewhere 



Finally, there are diagrams ^(i) and^(j), which contain C2 and C2 vertices. These vertices 
contribute for the first time at 0{kp) and introduce a dependence on the s-wave effective 
range and the p-wave scattering length a^. The diagrams are very similar to the one 
at 0{kp) and can be evaluated directly. One must be careful with the sign of the relative 
momenta k and k', since in the crossed contribution k' changes its sign. This is irrelevant 
for the C2 vertex, but for it changes the overall factor [g — 1) into {g + 1). Adding all 
contributions at 0{kp), we have 



£-^ = P77h (9-1) 



2M 



{kpasf kpTs + (^f + 1) {k^apf 



+ {g - 1){(0.07550 ± 0.00003) + - 3) (0.05741 ± 0.00002)} {kpa^ 



(26) 



which agrees with previous calculations |]T2|-|1J . 

Equation (|26|) illustrates that a conventional argument about the relative contribution 
of diagrams does not apply to the EFT-based diagrammatic expansion. Consider diagrams 
^(g) and |^(h) again. Diagram ^(g) has terms with two and four hole lines, while diagram 
^(h) has only terms with three hole lines and should be relatively suppressed according to the 

) shows that both diagrams are of the same order of magnitude, 



argument |^. Equation ( 
as expected from EFT power counting. We return to this point in the discussion below. 

In this section we have calculated the energy density of a dilute Fermi gas up to order 
{k-pR)^. In contrast to conventional calculations, the EFT approach is completely trans- 
parent and easily extended to higher orders, finite temperatures, and other observables. 
The EFT method is particularly advantageous when logarithmic divergences are present, as 
discussed in the next section. 



IV. LOGARITHMIC SINGULARITIES 

A naive extrapolation of Eq. (^) might lead one to conclude that the energy density 
has an analytic expansion in kp/A. This would be true if the Feynman graphs contained 
only power divergences. However, diagrams for three-to-three scattering have logarithmic 
divergences, and a consequence is that some terms in the energy density appear multiplied 
by powers of In/cp. Renormalization group methods provide a powerful tool to discern 
this analytic structure. We follow the discussion of Braaten and Nieto, who performed an 



analogous analysis for dilute Bose gases |]23 



As we discussed in Sect. ||, logarithmic divergences are special because they match 
onto logarithms of physical quantities. This means that we can identify where physical 
logarithms might appear by examining the renormalization group equations for the EFT 
coefficients. The general structure of these equations for dimensional regularization with 
minimal subtraction is 

= I3j{g) , (27) 
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FIG. 4. Diagrams for three-particle scattering in the vacuum that contribute to the running of 



where n is the (arbitrary) renormahzation scale, gj is a generic ra-body couphng [e.g., a C2i 
or a D2i), and Pj{g) is a polynomial in the couplings gi with A-independent coefficients. The 
terms in the polynomial correspond to diagrams which, if they have logarithmic divergences, 
will have nonzero coefficients. Matching A dimensions on each side of the equation [recall 
Eq. (pT])] greatly restricts the possible terms . 

The running of two-body couplings is determined by graphs in the two-body sector 
where there are only power divergences; therefore, none of these couplings depend on /i. So 
the lowest-dimension coupling that might run is Do, the coefficient of the three-body contact 
term in the EFT Lagrangian, 



1 



A-body = --Do{lp^1py 



(28) 



Since Dq ~ 1/MA^ and ~ 1/MA^*+^, Pj could have four Cq's or one Cq and one C2 or 
Cg. The latter two correspond to finite tree diagrams, so their coefficients in j3j are zero. 
That leaves only 2-loop diagrams with four Cq vertices, and an inspection of these uncovers 
logarithmic divergences in the diagrams in Fig. ^. 

Performing the loop integrals in D < 3 dimensions and then analytically continuing back 
to D = 3, the logarithmic divergences leads to simple poles l/{D — 3) in the 3 — 3 scattering 
amplitude. Braaten and Nieto have extracted the pole part of the diagrams in Fig. ^ 
(which involves removing an overlapping divergence). Using their result and taking into 
account the diagrams with cyclic permutations of the external momenta, the pole part of 
the 3 — s> 3 scattering amplitude is 



D-3 
1 



2(3-D) 



D~3 



- 2 In + 



(29) 



where the spin structure has been suppressed and the omitted terms vanish as D ^ 3. 
A counterterm removes the 1/{D — 3) pole in Eq. (^), leaving a finite term [Tj^ajfinite, 
which depends on the arbitrary scale /i. Since physical quantities cannot depend on /x, the 
renormalized coupling Dq must satisfy 



3— >3 1 finite 



iDf 



0, 



(30) 



or 
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where the hmit D ^ 3 has been taken.Q This has the form of Eq. (p^). 
Equation (^) is trivially integrated to find Dq^h): 

Doifi) = Do{l/as) + M' (Co)' ^^^^ ^''^'''^^ ' ^^^^ 

where we have chosen our reference scale to be the natural s-wave scattering length a^. 
Changes in the renormalization scale /i correspond to changes in how well the scattering 
at short distances is resolved. As less of the successive two-body scatterings in Fig. (|) is 
resolved, the contributions are shifted to the three-body contact term Dq. 

We can now use Eq. (|32D to obtain the contribution of order kp In kp to the energy 
density of the dilute Fermi gas. In order to obtain the full contribution at 0{kp), we 
would have to renormalize the 3 — 3 amplitude in the vacuum, perform the matching to 
determine Do{l/as) and then calculate all the energy diagrams at this order. (Alternatively, 
we could perform the matching with data from a many-body system.) To obtain just the 
nonanalytic contribution at this order, however, we have only to calculate diagram (k) in 
Fig. ^, since the logarithm in Eq. ( ^2]) must match to the In k-p from diagram (/) (plus others) 
26| . Equivalently, we keep only Fig. ^{k) with Dq evaluated at = kp, which avoids large 



cancellations between the diagrams. Applying the Feynman rules from the previous section 
we then have 

= p (^ - 2)(^7 - 1) H ^ (47r - 3v^) {kpa.r Hk^a,) , (33) 



which is the result obtained previously [p!^-[T^. Note that more than a two-fold degeneracy 



{g > 2) is needed to get a non-zero contribution, which follows from Pauli principle restric- 
tions. This procedure can be continued to higher orders to identify further (Infcp)" terms 
21. 



V. DISCUSSION AND SUMMARY 

We have derived the energy per particle in a dilute Fermi gas through terms of order 
kplnkp, reproducing results given in Refs. |T2|-|I3|. This calculation illustrates many of 
the virtues of an EFT approach to many-body physics. The EFT provides a model inde- 
pendent description of finite density observables in terms of parameters that can be fixed 
from scattering in the vacuum, with no off-shell ambiguities. Since the system is probed 
at low resolution, all vertices are contact interactions, which means exchange contributions 
are trivial and many diagrams vanish. After taking advantage of our freedom to choose a 
subtraction scheme, the power counting is exceptionally clean: each diagram contributes to 



Note that there are other diagrams contributing to ['Ta^alfinite- These diagrams, however, have 
no logarithmic divergences and consequently no /i dependence. 
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a single, well-defined order in the low-density expansion. Finally, the EFT provides tools for 
identifying and understanding the analytic structure of observables in terms of the short- 
(and long-) distance behavior. Future extensions of this work will include the generalization 
to finite temperature, single-particle properties, and linear response, as well as the complete 
renormalization at 0{kp) 

One can compare the clean, minimal calculation using an EFT to previous calculations, 
which required a more general set of diagrams that needed to be summed and expanded. 



with care taken to avoid double counting. In the standard treatments |]T2|-0, the focus 
is on the potential, rather than on observables as in the EFT. The EFT perspective says 
directly that at low energy the effective range parameters matter and not the details of a 
potential. A corollary is that one can use any complete set of regularized local interactions 
that reproduces the low-energy observables; there is no "best" potential. It's not that 
interactions are pointlike at low resolution, but that true observables are indistinguishable 
(up to a well-defined truncation error) from those calculated with pointlike interactions. 
These interactions are wrong at short distances, but renormalization theory prescribes how 
the discrepancy from incorrect high energy behavior can be systematically corrected by the 
renormalization of the EFT coefficients. 

The idea of replacing the "true" potential with one that sums low-energy scattering effects 
and is easier to work with (a "pseudo-potential") is an old idea |2^. The EFT framework 



makes it systematic and extendable. Even if we did know the underlying interaction {e.g., 
a hard sphere potential), it is simpler and more efficient to match in the vacuum to an EFT 
and then to use the EFT for the finite density calculation. This is analogous to the use 
of effective theories for heavy-quark physics (NRQCD) and QED bound states (NRQED), 
which are beneficial even though the underlying theories are known (cf. Ref. |l|]). 

The EFT analysis can also clarify the nature of many-body systems. In Ref. the 



energy expansion for a dilute Fermi gas is followed by a discussion of the hole-line expansion, 
which is used to do a form of power counting. The simple argument is that each additional 
hole-line is accompanied by momenta restricted by fcp, leading to an additional factor of 
density that suppresses the diagram at low density. Particle lines have unrestricted momenta 
and are therefore not suppressed, and so must be included to all orders. The power counting 
in our EFT reveals that a hole-line expansion for the low-density natural Fermi gas is only 
an artifact of a subtraction scheme that does not remove power divergences. Indeed, the 
particle-particle and particle-hole ring diagrams in Fig. ^ contribute to the same order in kp 
and are numerically comparable. This is cleanly reproduced in the EFT because of minimal 
subtraction, which results in no distinction between particle and hole lines in the power 
counting. 

We saw that going further in the dilute Fermi gas energy expansion requires input from 
three-to-three scattering (or other many-body observables). This means that three-body 
(and four-body) contributions are inevitable. If only two-body potentials were used, one 
would find different results at this order in the expansion for different potentials that agreed 
for two-body scattering, which would conventionally be attributed to a lack of knowledge of 
the off-shell two-body potential. From the EFT perspective, the differences would be due 
to short- distance behavior that must be corrected with many-body counterterms. For the 
three-nucleon system, accounting for the three-body contact interaction leads to a compelling 
explanation of the Phillips line The need to renormalize because of incorrect high- 
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energy behavior is not avoided because of finite loop integrals {e.g., due to a form factor). 
A sensitivity to cutoff parameters is a signature that something is missing in the theory. 

The extension of the EFT procedure for dilute, repulsive Fermi systems to the nuclear 
matter problem [^,0 is far from immediate. In particular, the EFT subtraction and 
renormalization scheme used here is not appropriate for a treatment of nuclear matter. 
That is because the EFT momentum breakdown scale will be set by the lowest mass scale in 
the effective range expansion, which is of order 10 MeV because of the fine-tuned scattering 
lengths. Since nuclear matter equilibrium density is around kp ^ 270 MeV, the expansion 
parameter kp/A ~ 27 and the perturbation series is not useful. Alternative power counting 



schemes have been proposed to deal with this problem in the vacuum |10,11|. The consistent 



extension of such a scheme to finite density, which may reintroduce a hole-line expansion, is 
the next step. However, an EFT with only short-range interactions will still fall far short, 
with a breakdown scale of 771,^/2 [Q. Thus, pions will have to be included as long-range 
physics. We are confident that the virtues of the EFT methods illustrated for a much 
simpler system here will continue to pay dividends as they are adapted to the study of 
nuclear matter. 
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